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ABSTRACT 


A  theoretical  study  has  been  made  of  the  scattering  of  electro 
magnetic  waves  from  both  underdense  and  overdense  turbulent  low  alt 
tude  rocket  plumes.  The  analysis  considers  a  number  of  parametric 
assumptions  about  the  nature  of  the  turbulence,  namely ,  type  of  cor¬ 
relation  function,  correlation  lengths,  isotropy  of  turbulence, 
finite  conductivity  and  plume  intermittency .  The  analytical  re¬ 
sults  have  been  applied  to  the  cross  sections  of  both  an  underdense 
and  overdense  plume  at  65,000  ft  altitude.  The  most  important  re¬ 
sult  is  that  the  theoretical  cross  section  can  be  large  for  some 
reasonable  assumptions  about  the  nature  of  the  plume  turbulence. 


INTRODUCTION 


Analytical  and  experimental  study  of  the  effects  of  radar  re¬ 
flection  from  plumes  has  been  an  area  of  active  research  for  many 
years.  In  this  paper  we  concentrate  on  the  analytical  aspects  of 
the  problem.  Our  plan  is  to  formulate  models  that  bear  a  reason¬ 
able  resemblance  to  a  practical  plume  scattering  problem  and  yet 
are  simple  enough  to  handle.  We  will  introduce  into  the  models 
some  scattering  phenomena  that  have  not  been  introduced  into  prac¬ 
tical  rocket  plume-EM  wave  scattering  problems  before. 

We  do  not  mean  to  say  that  these  are  the  only  phenomena  that 
should  be  considered  in  the  electromagnetics  of  plume  scattering. 
They  are  not.  In  fact,  the  numerical  results  of  this  paper  (i.e., 
the  calculated  cross  sections)  vary  over  a  wide  range  of  values 
depending  upon  the  nature  of  the  model  of  the  scattering  process. 
However,  if  we  are  to  make  progress  in  this  very  complicated  field 
we  must  choose  a  large  number  of  models  for  our  theoretical  analy¬ 
sis  to  compare  with  experimental  data.  By  a  long  and  tedius  process 
one  will  arrive  at  models  that  correctly  describe  the  electromag¬ 
netics  under  all  operating  conditions. 

The  analytical  studies  of  the  problem  have  in  general  been 
concerned  with  two  main  areas  of  interest,  namely  fluid  mechanics 
and  electromagnetics.  The  fluid  mechanics  approach  considers 
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the  problem  of  determining  the  electron  density,  density  fluctua¬ 
tion,  temperature  and  collision  frequency  distribution  throughout 
the  plume.  Many  analytical  models  have  been  concerned  with  the 
various  gas  dynamic  effects  that  could  influence  the  electron 
properties,  some  of  which  are  both  equilibrium  and  nonequilibrium 
chemistry,  turbulence  in  the  plume  core,  viscous  interaction  ef¬ 
fects,  diffusion  of  the  plume  into  a  moving  environment,  tempera¬ 
ture  differences,  and  inhomogeneity. 

The  study  of  the  electromagnetic  effect  has  generally  been 
separated  into  two  categories.  The  first  is  the  volume  scattering 
problems  where  the  illuminating  energy  to  first  order  passes  un¬ 
attenuated  through  the  plume.  The  scattered  energy  results  from 
the  electrons  excited  by  this  energy.  Secondly,  the  problems  can 
be  considered  as  surface  scattering  problems  where  the  plume  elec¬ 
tron  concentration  is  high  enough  that  one  can  define  a  surface 
like  a  metal  plate  from  which  scattering  will  take  place. 

In  this  paper  we  will  concentrate  on  the  study  of  the  electro¬ 
magnetic  aspects  of  the  problem.  We  will  briefly  describe  the 
analysis  that  was  made  on  the  fluid  mechanics  aspects  and  show 
how  the  results  of  this  turbulent  calculation  were  used  as  a  "take 
off"  point  for  the  description  of  the  analytical  models  of  the 
turbulent  plume  used  in  the  EM  wave  scattering  analysis. 
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PLUME  CALCULATION  FOR  VOLUME  SCATTERING  PROBLEM 


Two  computer  programs  were  used  to  determine  the  electron 
distribution  in  the  plume.  The  first,  a  nozzle  program,  cal¬ 
culated  the  ionized  and  neutral  species  concentration  in  the  com¬ 
bustion  chamber  and  the  nozzle.  The  most  important  ionizing 
species  of  the  liquid  propellant  are  the  alkali  metal  impurties 
of  sodium  and  potassium.  Their  molar  concentration  is  about 
10  ppm.  A  second  computer  program  is  used  for  the  fluid  analysis 
of  the  rocket  exhaust  plume.  The  species  concentration  at  the 
nozzle  exit  plane,  as  calculated  by  the  previous  program,  is  used 
as  initial  condition  for  this  program.  The  program  considers 
turbulent  mixing  and  nonequilibrium  chemistry  in  the  calculation 
of  the  neutral  and  ionized  species  concentrations  throughout  the 
plume  that  is  flowing  into  a  moving  environment.  The  plume  pro¬ 
gram  is  a  Grumman  adaptation  of  an  existing  program  (Ref.  1). 

The  distribution  of  the  electron  concentration  obtained  from  the 
exhaust  plume  computer  calculation  is  used  to  construct  an  analyti¬ 
cal  model  of  the  plume  that  has  a  relatively  simple  geometric  shape. 

We  have  sketched  the  actual  electron  concentration  of  a  liquid 
propellant  axisymmetric  plume  from  the  computer  print-out  of  the  ex¬ 
haust  plume  program  (Fig.  1).  These  results  are  for  an  altitude  of 
65,000  feet.  Although  not  shown  on  the  figure,  the  overdense  region 
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of  the  plume  for  the  frequencies  considered  is  very  small;  hence 
scattering  from  this  region  will  be  neglected.  On  the  other  hand, 
there  appears  to  be  a  very  large  region  of  weakly  ionized  gas  ex¬ 
tending  as  far  as  150  meters  from  the  nozzle  exit,  with  a  maxi¬ 
mum  radius  of  about  10  nozzle  exit  radii.  From  this  observation 
the  first  conclusion  that  we  can  draw  about  this  liquid  propellant 
plume  is  that  it  will  be  underdense.  Therefore  we  must  consider  a 

volume  scattering  analysis.  Further  investigation  shows  that  the 

12 

electron  density  varies  along  the  axis  from  approximately  10 

at  the  nozzle  exit  plane,  decreasing  very  rapidly  at  first,  and 
9 

reaching  10  at  about  150  meters  from  the  nozzle  exit  plane. 

The  radial  variation  at  every  axial  position  varies  very  slowly 

from  a  maximum  value  on  the  axis  out  to  some  radial  position  and 

then  drops  sharply  to  a  relatively  low  value.  The  curve  has  a 

knee  shape.  As  an  example  of  this,  consider  the  radial  variation 

at  an  axial  location  of  28  meters.  The  electron  density  is 
10  9 

1 . 3  x  10  on  the  centerline,  10  at  five  nozzle  radii  and 
drops  sharply  to  1.8  x  10^  at  6.25  nozzle  radii.  The  radial 
position  of  the  knee  varies  in  an  approximately  conical  manner 
from  the  nozzle  exit  plane  to  an  axial  position  about  28  meters 
from  the  nozzle  exit  plane.  The  position  of  the  knee  remains 


relatively  constant  at  about  10  nozzle  radii  from  that  point  to 
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150  meters  from  the  nozzle  exit  plane.  One  could  consider  the 

electron  distribution  as  contained  in  a  cylinder  of  radius  10 

i.ozzle  exit  radii  and  150  meters  long,  with  the  electron  den- 

12  9 

sity  decreasing  monotonically  from  10  to  10  along  the  axis. 

For  the  volume  scatter jog  problem  we  have  taken  the  model  of  the 
turbulent  plume  to  be  a  cylinder  150  meters  by  7.5  nozzle  exit 
radii  in  radius  with  an  average  uniform  electron  distribution. 

The  axial  variation  of  the  electron  concentration  has  been  neglected. 

ELECTROMAGNETICS  OF  VOLUME  SCATTERING 

Volume  scattering  of  electromagnetic  waves  from  turbulent  fluid 
plasmas  was  first  considered  by  Booker  and  his  associates  (Ref.  2) 
in  the  early  Fifties  and  has  been  considered  by  others  in  a  similar 
manner  since  that  time  (Refs.  3  and  4).  Physically,  Booker's  assump¬ 
tion  is  that  the  plasma  is  essentially  transparent  to  the  primary 
radiation  from  the  external  source  .  The  scattered  or  reradiated 
energy  is  generated  by  the  excitation  of  the  electrons  in  the  turbu¬ 
lent  plasma  by  this  primary  radiation.  The  mathematical  technique 
describing  this  scattering  mechanism  is  called  the  Born  approxima¬ 
tion  and  was  first  developed  to  describe  quantum  mechanical  scat¬ 
tering  phenomena.  Under  the  assumption  of  the  Born  approximation 
(where  multiple  scattering  is  not  important)  we  can  write  d  vn  the 
integral  for  the  scattered  power,  which  is  essentially  nothing  more 
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than  a  solution  to  Maxwell's  equations  for  the  scattered  radiative 
power  for  a  collection  of  sources  (excited  electrons)  in  some 
volume  in  space.  The  cross  section,  o,  represents  the  ratio  of 
the  scattered  power  in  a  certain  direction  to  the  incident  power. 
This  integral  is 


a 


2  2 
sin  a  7 

e 


(1) 


-(ki-kr)  . 


dV(1)dV(2) 


The  time  average  of  the  electron  density  at  two  points  appears 
within  the  integral 


k  ,k 
i*  r 


-  incident,  reflected  wave  vectors 


classical  electron  radius 


V(D,V(2) 


angle  between  the  incident  electric  field 
vector  and  the  scattered  wave  vector 

volume;  superscripts  represent  two  inde¬ 
pendent  integrations 


The  term  <  N2  )  can  be  expanded  into  the  form 
<  Nx  N2  )  -  Nl0  N20  +  <  ANX  £N2  ) 


(2) 
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where  ^q,N2q  are  the  time  average  or  mean  values  of  the  elec¬ 
tron  concentration  at  the  point  1,2,  and  ^^>^2  are  the 
fluctuating  values  of  the  electron  concentrations.  When  Eq.  (2) 
is  substituted  into  Eq.  (1),  the  scattering  cross  section  is  com¬ 
posed  of  two  terms:  the  first  is  the  scattering  from  the  mean 
electron  concentration,  called  the  coherent  scattered  energy;  the 
second  is  the  scattering  from  the  fluctuating  part,  called  the 
incoherent  scattered  energy.  We  find  that  in  most  rocket  plume-F.M 
wave  volume  scattering  problems  where  the  fluid  flow  is  turbulent, 
the  incoherent  scattered  energy  is  by  far  the  greater. 

The  usual  procedure  in  calculating  the  volume  scattering 
cross  section  of  the  turbulent  body  of  ionizeJ  gas  is  to  intro¬ 
duce  the  two  point  autocorrelation  functions  of  the  turbulence 
for  the  term  <  AN-^  AN2  )  which  can  be  integrated,  either  nu¬ 
merically  or  analytically.  The  form  of  the  autocorrelation  func¬ 
tion  is  such  that  it  decreases  monotonically  as  the  distance  in¬ 
creases  between  the  spatial  points  of  AN^  and  AN2  .  The  char¬ 
acteristic  length  of  this  decay  is  called  the  correlation  length. 
The  autocorrelation  function  has  the  form 


(  AN1(x1)  AN2(x2)  >  «=  <  (INj^)2  )f(x1-x2)  .  (3) 

An  example  is  the  gaussian 

(xj^  -  x2)2 
'  .  2 

c  t  \  2a 

f(xx-  x2)  =  e 

where  a  is  the  correlation  length.  In  the  volume  scattering 

problem  treated  here  we  will  consider  only  isotropic  turbulence. 
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Therefore  the  function  f  will  vary  only  as  the  magnitude  of  the 
distance  between  the  spatial  points. 

The  usual  procedure  in  volume  scattering  problems  is  to  assume 
that  the  characteristic  length  of  the  turbulence  is  much  smaller 
than  the  smallest  linear  dimension  of  the  plasma.  If  one  considers 
this  case  and  substitutes  Eq .  (3)  into  Eq.  (1)  (considering  the  in¬ 
coherent  scattering  only)  and  integrates  over  the  volume  of  the 
ionized  gas  one  finds 

?  =  (sin2C)(^)Vi(f,a)  .  (4) 

In  Eq.  (4)  we  have  substituted  a  discrete  function  for  the 

autocorrelation  function  f  (for  example,  the  gaussian  or  the 

exponential).  Since  the  correlation  function  is  expressed  in  terms 

of  a  difference  between  two  points  in  a  gas,  a  coordinate  trans- 

1  .  1 

formation  has  been  made  substituting  sum  (X2)  and  dilference  x^ 
for  the  real  coordinates 


1 

x^  =  x^  -  X2 

(5) 

x2  =  ^X1  +  x2  ^ 

Because  the  correlation  length  is  small  relative  to  the  plasma 
size,  the  function  f  dies  out  very  fast  for  large  values  of  the 
di i ference  coordinates.  Therefore  the  limits  of  integration  for 


S 


This  integral 


the  difference  coordinates  can  be  replaced  by  ±*. 
then  simply  becomes  the  three  dimensional  spectral  density  £(f,a) 
of  the  particular  autocorrelation  function  f  with  its  particular 
correlation  length  a.  Equation  (4),  sometimes  called  the  Booker- 
Gordon  relation,  is  the  conventional  cross  section  relation  to  use 
when  volume  scattering  in  the  Born  approximation  is  applicable  and 
when  the  scale  of  the  turbulence  is  small  compared  to  the  size  of 
the  plasma. 

However,  recent  experimental  measurements  indicate  that  the 
turbulence  scale  may  not  always  be  small  compared  to  the  size  cf 
the  plume.  Slattery  (Refs.  5  and  6)  and  his  associates  at  Lincoln 
Laboratory  have  made  microdensitometer  tracings  from  schlieren 
photographs  of  light  passing  through  the  hypersonic  turbulent  wake 
of  spheres  fired  in  a  ballistic  range.  Using  a  mathematical  tech¬ 
nique  developed  by  Uberoi  and  Kovasznay,  they  were  able  to  reduce 
the  correlation  in  the  contrast  on  the  photographic  plate  to  an 
actual  autocorrelation  function  of  the  turbulence  in  the  wake. 

These  results  show  that  the  correlation  length  (defined  as  the  dis¬ 
tance  required  for  the  correlation  function  to  reduce  in  value  1/e) 
is  constant  and  approximately  one-half  the  diameter  of  the  sphere 
generating  the  wake.  Figure  2  is  taken  from  these  results  and 
shows  that  the  mean  values  of  the  wake  diameter  to  correlation 
length  ratio  can  vary  from  1  at  10  sphere  diameters  behind  the 
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wake  to  10  at  1000  sphere  diameters  behind  the  wake  for  a 
3/8  in.  sphere. 

Measurements  of  statistical  turbulence  properties  have  been 
made  at  Stanford  Research  Institute  (Ref.  7)  by  seeding  an  ethylene 
oxygen  flame  with  potassium  chloride.  In  this  case,  local  measure¬ 
ments  were  made  of  the  fluctuation  cf  the  ion  and  electron  density 
with  ion  and  electron  probes.  Correlation  coefficients  were  de¬ 
termined  from  the  probe  measurements  by  a  signal  correlator.  The 
results  of  their  data  taken  from  the  report  are  shown  in  Fig.  3. 
Mean  values  of  the  plume  diameter  to  correlation  length  ratio 
are  4.0  at  9.0  nozzle  diameters  downstream  and  5.3  at  20.0 
nozzle  diameters  downstream  of  the  nozzle  exit  plane.  The  above 
two  experiments  indicate  that  under  some  physical  conditions  the 
turbulence  scale  lengths  of  wake  and  jets  may  not  be  small  compared 
to  the  size  of  the  exhaust  plume.  This  must  be  considered  in  our 
electromagnetic  scattering  analysis.  A  full  discussion  of  the 
analysis  for  this  case  is  presented  in  Ref.  8.  We  will  sketch  only 
the  important  highlights  here,  showing  that  for  a  large  correla¬ 
tion  length  the  finite  volume  cross  sections  can  be  large. 

The  finite  volume  effect  is  easily  seen  in  a  one  dimensional 
problem.  Substituting  Eq.  (3)  in  the  incoherent  scattering  part 
of  Eq.  (1)  in  one  dimension  only  (Fig.  4),  we  find  that  the  colli¬ 
sion  cross  section  can  be  written  in  the  form 
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a  *  2A 


n2A  1 

1 

1  ~ikxl 
g(xp  e 

_  O  A 

1  2A 

k.  -  k  | 

I  r 1 


dx. 


(6) 


1  2  1 

where  g(x^)  »  <  (AN)  )  f(x^),  the  turbulence  correlation  func- 

2  2  1 

tion  and  suppressing  sin  n  The  x^  is  the  difference  co¬ 

ordinate  [as  in  Eq.  (5)].  Use  of  sum  and  difference  coordinates 
amounts  to  a  rotation  of  the  coordinate  system  (Fig.  5).  Integra¬ 
tion  will  proceed  from  line  4  to  line  1  and  from  line  3  to  line  2 


1 

;2 


in  the  xX  coordinate.  The  lines  are  then  summed  in  the  x 


difference  coordinate.  If  we  let  x^  -  ay,  Eq .  (6)  becomes 

p2A/a 


a  =  2A 


ag(y)  e 


■ikay 


1  - 


2A/a 


dy 


(7) 


-2A/a 


If  we  now  let  A/a  -*  «,  the  second  term  in  the  bracket  approaches 
zero  (A/a  is  in  the  denominator)  and  the  first  integral  takes 
the  form 


a 


2A 


,  1.  -ikx  ] 
g(x  )  e  dx 


(8) 
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Equation  (8)  is  the  one  dimensional  Booker-Gordon  equation,  the 
integral  being  the  spectral  density  of  the  autocorrelation  func¬ 
tion  and  2  A  being  the  one  dimensional  "volume."  However,  when 
A/a  is  finite  we  see  that  an  extra  term  appears  in  the  scattering 
integral  [Eq.  (7)].  Under  these  conditions,  incorrect  results 
would  be  obtained  if  the  Booker-Gordon  integral  were 

The  problem  of  evaluating  the  scattering  cross  section  for  a 
three  dimensional  finite  volume  plasma  is  a  generalization  of  the 
evaluation  of  the  scattering  cross  section  for  the  one  dimensional 
plasma.  A  set  of  transformation  equations  is  introduced  [similar 
to  Eq.  (5)1  for  each  of  the  three  real  space  coordinates.  Since 
we  have  to  integrate  over  the  volume  twice,  as  in  the  one  dimen¬ 
sional  problem,  the  integration  will  be  over  a  six  dimensional 
volume.  The  scattering  cross  section  will  be  represented  by 


a  = 


•  2o  2 

sin  P.  y 

e 


g(x}» 


1  K 

yl’Zl> 


ik 


dx.  dy,  dz 


1 


(9) 


The  subscripts  1,2  are  the  difference  and  sum  coordinates,  re¬ 
spectively,  as  in  Eq .  (5).  The  limits  of  integration  will  con¬ 
tain  the  surface  of  the  plasma  as  in  the  one  dimensional  scattering 
relation.  Since  we  are  considering  a  homogeneous  turbulence,  the 
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correlation  length  will  be  a  function  of  Che  difference  coordinates 
only.  If  we  consider  the  integration  in  our  one  dimensional  plasma, 
again  we  see  that  it  is  really  in  two  regions  of  the  transformed 
space  (Fig.  5).  The  limits  of  integration  for  the  double  integral 
are  different  in  each  region.  Therefore  the  one  dimensional  scat¬ 
tering  cross  section  is  really  a  sum  of  two  integrals  in  transformed 
space.  For  the  three  dimensional  plasma  the  scattering  cross  sec¬ 
tion  will  be  a  sum  of  eight  integrals  over  eight  regions  in  the 

six  dimensional  space.  We  can  write  this  sum  of  eight  integrals  as 

i,j,k-2 

a  =  sin2fi  %2  ^  I,  .  .  (10) 

e  l.  k  j  1  v  ' 

i > j ,k-l 


One  of  the  integrals  can  be  written  as 
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1 
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J  c 

c 
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t 

t 

.,1,1.  1,1,  1,1 

h  dx£  dx^  dy2  dy^  dz^  dz^ 


(U) 


4  3  4 

where  the  upper  limits  are  the  subscripts  of  l^j  i  (1^)  • 

The  numbers  over  the  integrals  represents  the  limits  of  integra¬ 
tion  from,  say,  line  4  to  line  1  in  the  x  and  z  coordinates 
and  line  3  to  line  2  in  the  y  coordinate.  The  integration  pro¬ 
ceeds  in  the  order  x,y,z.  In  the  one  dimensional  case  the  equa¬ 
tion  tor  lines  1  through  4  contains  the  parameter  ±A,  the  boundary 
of  the  plasma.  In  the  case  of  the  three  dimensional  plasma,  lines  1 
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through  4  also  contain  the  boundary  of  the  plasma;  however,  they 
are  not;  in  a  functional  form.  They  are  called  where 

/,tn,n  vary  from  1  to  2  and  are  known  functions  of  the  geometry 
of  the  plasma  (Ref.  8).  The  equations  for  the  limits  of  integra¬ 
tion  are  numbered  in  a  similar  manner  as  for  x  as  shown  in  Fig.  5. 
Therefore  we  see  that  there  are  two  integration  regions  for  each 
coordinate,  which  make  a  total  of  eight  combinations  of  six  dimen¬ 
sional  integrations;  thus  the  eight  regions  in  the  six  dimensional 
space . 

The  actual  limits  of  integration  are  given  in  Ref.  8,  and  a 
table  of  their  values  is  reproduced  in  Fig.  6.  We  might  note  here 
that  the  eight  integrals  will  contract  to  the  single  Booker-Gordon 
scattering  integral  when  the  correlation  length  is  small  compared 
to  the  smallest  characteristic  dimension  of  the  plasma.  The  above 
relation  holds  for  any  shaped  three  dimensional  plasma. 

The  Booker-Gordon  analysis  shows  that  the  cross  section  is 
independent  of  the  shape  and  the  orientation  of  the  plasma  rela¬ 
tive  to  the  transmitter  and  receiver,  but  depends  only  on  the 
plasma  volume.  The  finite  volume  analysis  shows  that  the  orienta¬ 
tion  and  the  shape  of  the  plasma  influence  the  value  of  the  scat¬ 
tering  cross  section. 

In  theory  it  is  possible  to  calculate  the  scattering  cross 
section  of  any  homogeneous  plasma  shape,  with  any  correlation 
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function  by  means  of  Eq.  (10)  with  the  limits  of  integration  given 
in  Fig.  6.  However,  the  evaluation  of  the  eight  integrals  for  a 
correlation  function  that  is  not  separable  could  be  a  lengthy  pro¬ 
cess.  The  first  finite  volume  that  was  considered  was  a  parallel¬ 
epiped  plasma  (Fig.  7)  using  a  gaussian  correlation  function. 

Under  these  conditions  the  are  constants  (the  boundaries 

of  the  plasma,  i.e.,  -A,  A,  -B,  B,  -C,  C) .  The  gaussian  correla¬ 

tion  function  is  separable  in  the  x,y,  and  z  coordinates.  Since 
the  limits  of  integration  for  the  boundaries  are  constants  and  the 
correlation  function  itself  is  separable,  it  i':  possible  to  perform 
the  integration  first  in  the  x  sum  and  difference  coordinate, 
then  in  the  y  sum  and  difference  coordinate,  and  finally  in  the 
z  sum  and  difference  coordinate.  Then  each  of  the  sum  coordinates 
can  be  integrated  so  that  the  cross  section  is  expressed  as  a  prod¬ 
uct  of  three  single  integrals  in  the  x,y,  and  z  difference  co¬ 
ordinate.  We  can  write  the  cross  section  as 

a  =  KxK2K3  (12) 

where 


There  are  similar  expressions  for  K2  and  with  the  limits 

of  integration  being  -2b,  2B  and  -2C,  2C ,  respectively.  We  can 


15 


see  that  the  cross  section  for  the  rectangular  parallelepiped  is 
just  a  product  of  three  one  dimensional  cross  sections.  We  can 
understand  the  physics  of  the  scattering  mechanism  by  studying 
just  >'-?  one  dimensional  cross  section.  If  we  consider  the 
gaussian  correlation  function 


and  make  the  substitution 


‘v  2  a 


the  relation  for  becomes 


3/2 

aFy  -  2  /<£Aa 


£ 


0 


(1  -  p  e-y  cos  V?  kay  dy 


(13) 


where  3  =  ^  2A /a  and  the  one  dimensional  Booker-Gordon  relation  is 


2,  2 

a  k 

0  3  /  2  i  '  2 
agG  «=  2  Aav2  e 

The  ratio  of  finite  volume  to  Booker-Gordon  cross  section  is 

.2  3 
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Equation  (15)  has  been  solved  numerically  on  a  computer  and  the 
results  have  been  plotted  (Figs.  8  and  9).  It  can  be  seen  from 
the  figures  that  for  ka  >  3,  the  cross  section  can  be  very  large 
compared  to  the  Booker-Gordon  values  for  small  values  of  A/a. 

VOLUME  SCATTERING  FROM  A  PLUME  AT  65,000  FEET 


The  fluid  mechanic  and  electromagnetic  techniques  outlined 
in  the  previous  section  have  been  applied  to  the  problem  of  deter¬ 
mining  the  backscatter  scattering  cross  section  of  a  turbulent 
rocket  exhaust  plume  at  65,000  ft  altitude.  The  turbulent  plume 
model  considered  was  150  meters  long  with  a  radius  of  7.5  noz¬ 
zle  radii.  The  turbulence  was  considered  to  be  isotropic  and  homo¬ 
geneous  with  a  gaussian  correlation  function.  For  the  finite  vol¬ 
ume  correction  to  the  scattering  cross  section  we  use  a  rectangu¬ 
lar  parallelepiped  with  a  side  dimension  of  11  meters  (2A  «  11m). 
Calculations  were  made  with  three  assumed  values  of  the  plume  cor¬ 
relation  lengths  of  0.0055,  0.055,  and  0.55  m.  The  correspond¬ 
ing  values  of  A/a  are  then  1000,  100,  and  10.  Cross  sections 
were  calculated  with  three  orientations  of  the  plume  relative  to  the 
radar.  The  orientations  were  normal  to  the  plume  axis,  along  the 
plume  axis,  and  at  45°  to  the  plume  axis.  Four  radar  frequencies 
were  considered:  200  MHz,  400  MHz,  3GHz ,  and  10GHz .  The  results 
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of  che  calculations  are  given  in  Fig.  10.  Where  finite  volume 
corrections  were  made,  we  see  this  indicated  on  the  figure  by  an 
asterisk.  Finite  volume  corrections  are  needed  when  the  value  of 
ka  exceeds  1.  Normalized  scattering  cross  sections  vary  from  a 
high  value  of  2660  to  a  low  value  of  0.03,  depending  on  the 
parameters  considered.  The  cross  section  can  be  large  for  some 
parametric  assumptions  about  the  nature  of  the  plume  turbulence. 

PLUME  CALCULATION  AND  ELECTROMAGNETICS 
FOR  A  SURFACE  SCATTERING  PROBLEM 

The  same  two  computer  programs  were  used  to  calculate  the 
electron  distribution  in  the  plume  for  the  surface  scattering  prob¬ 
lem  as  for  the  volume  scattering  problem.  The  reaction  and  rate 
constants  used,  however,  were  for  the  solid  propellant.  The  nozzle 
program  also  predicted  the  formation  of  some  aluminum  oxide  in  the 
exhaust  products.  This  material  will  form  a  particulate  matter, 

95  percent  of  which  will  be  less  than  10  microns  diameters  in 
size.  As  these  particles  flow  from  the  combustion  chamber  and  out 
the  nozzle,  their  inertia  prevents  them  from  spreading  radially 
outward  to  the  edge  of  the  nozzle.  The  heavier  particles  remain 
close  to  the  axis  of  the  nozzle.  The  study  also  shows  that  the 
particle  density  at  the  exit  plane  is  relatively  sparse,  varying 
from  a  high  value  of  2.2  particles/cc  on  the  centerline  for  the 
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1.5  micron  size  to  about  0.018  particles/cc  at  a  distance  of 
about  half  an  exit  plane  radius  for  the  8-micron  size.  One 
would  expect  these  low  density  and  particle  inertia  effects  to 
continue  out  the  nozzle.  Therefore  a  relatively  narrow  cone  of 
low  density  particulate  material  would  be  seen  extending  out  into 
the  plume.  Since  the  particulate  matter  is  a  dielectric  and  its 
number  density  is  very  low  throughout  the  plume,  the  electromag¬ 
netic  scattering  cross  section  of  this  material  will  be  very  low. 
Therefore,  electromagnetic  scattering  from  the  particulate  matter 
is  neglected  in  this  investigation. 

The  alkali  metal  concentration  in  the  solid  propellant  was 
about  200  ppm  and  this  made  the  electron  concentration  in  the 
plume  about  20  times  greater  than  that  of  the  liquid  propellant 
calculation.  We  are  now  interested  in  the  scattering  cross  sec¬ 
tion  for  400  MHz  waves.  These  two  facts  combined  to  make  a  large 
axisymroetric  overdense  region  in  the  plume.  This  volume  expands 
rapidly  to  4.6  meters  in  diameter  close  to  the  nozzle  exit  plane, 
then  bows  out  gradually  to  6.5  meters  and  finally  tapers  to  2.6 
meters  in  diameter  at  150  meters  from  the  nozzle  exit  plane.  The 
average  plasma  frequency  in  the  overdense  volume  is  10^  radians/ 
second  and  the  radar  frequency  is  2 . 5  x  10  radians/second.  Since 
the  overdense  volume  is  so  large  for  the  solid  propellant  rocket  at 
this  altitude,  we  assume  that  the  reflection  from  this  fluctuating 
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surface  is  the  dominant  process  in  this  analysis.  A  sketch  of  the 
overdense  volume  as  actually  computed  by  the  exhaust  plume  program 
is  shown  in  Fig.  11. 

Based  on  the  calculated  overdense  volume  described  above,  the 
following  analytical  model  for  the  surface  scattering  is  proposed. 
The  overdense  volume  is  approximated  as  a  cylinder  of  150  meters 
in  length  and  4.6  meters  in  diameter.  Corrections  to  the  scat¬ 
tering  cross  section  are  then  made  for  a  finite  conductivity 
plasma  and  for  surface  roughness  due  to  turbulence.  These  correc¬ 
tions  are  made  by  multiplying  the  smooth  infinite  conducting  sur¬ 
face  cross  sections  by  correction  factors  for  finite  conductivity 
and  surface  roughness.  Thus  the  corrected  expression  for  the 
cross  section  would  be 


'  <cf=>  •  <V 


(16) 


where  oc  is  the  corrected  cross  section,  agc  the  smooth  cylinder 
cross  section,  C^c  the  finite  conductivity  correction  and  CR  the 
roughness  correction.  The  analytical  expression  used  is  for  the 
backscatter  cross  section  at  broadside  for  a  cylinder  where  the 
length  and  radius  are  large  compared  to  the  radar  wavelength.  The 
cross  section  expression  for  a  smooth  cylinder  is  given  by  the 
following  expression  (Ref.  9) : 

2 


a  =  4  k  Rt 
sc  o 


(17) 
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where  R  is  the  cylinder  radius  and  t  is  the  half  length  of 
the  cylinder. 

The  finite  conductivity  correction  is  based  on  a  planar 
plasma-free  space  interface.  Finite  conductivity  arises  as  a 
result  of  collisions  of  electrons  with  neutral  particles  in  the 
plume.  The  average  collision  frequency  in  the  overdense  region 
of  the  plume  is  1.2  x  10^.  One  method  for  correction  of  the 
finite  conductivity  effects  is  proposed  by  Hochstim  (Ref.  10). 

He  calculates  the  reflection  coefficient  of  a  semi-infinite 
ionized  gas  medium  with  a  plane  boundary.  Figure  12  is  taken 
from  Hochstim's  report.  It  can  be  seen  that,  tor  the  conditions 
of  this  problem  (vc“  1.2  x  10^,  -  10^,  to  «  2.5  x  10^), 

the  reflection  coefficient  is  about  0.32.  Since  the  cross  sec¬ 
tion  correction  is  proportional  to  the  square  of  the  reflection 
coefficient,  the  finite  conductivity  correction  factor  (C^) 
would  be  0.1. 

Hochstim's  analysis  is  for  a  semi-infinite  medium.  French, 
Cloutier,  and  Bachynski  (Ref.  11)  have  calculated  the  reflection 
coefficient  for  the  case  of  a  slab  2.5  plasma  wavelengths  thick 
and  also  for  a  semi-infinite  slab.  For  the  conditions  of  our 


plasma  (v  /a>  **  1)  there  is  little  difference  between  the  re- 

c  P 

flection  coefficients  for  the  two  thicknesses  given  above.  Only 


for  small  values  of  v  /o>  (less  than  0.1)  where  there  are  wide 

c  p 
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periodic  variations  in  the  reflection  coefficient  as  oo/co  in¬ 
creases  for  the  plasma  of  thickness  2.5  X  would  it  be  possible 
to  have  large  differences  in  reflection  coefficient  for  the  thick 
and  thin  slab  cases.  For  the  average  plasma  and  collision  fre¬ 


quencies  of  the  overdense  volume,  French  et  al.  determine  a  re¬ 


flection  coefficient  of  0.3  or  a  correction  to  the  cross  section 


of  0.09.  Since  the  thin  slab  and  semi-infinite  volume  finite 


conductivity  corrections  are  almost  the  same  value,  one  could 
assume  that  the  correction  factor  for  the  cylinder  is  0.1. 

The  second  correction  factor  to  be  applied  to  the  cross  sec¬ 
tion  for  the  smooth  cylinder  is  that  of  surface  roughness.  This 
roughness  arises  because  of  the  turbulent  nature  of  the  flow  and 
manifests  itself  as  a  moving,  or  undulating,  surface  where  the 
motion  of  neighboring  positions  is  statistical  in  nature.  One 
analytical  approach  that  one  could  take  is  to  carry  over  to  our 
problems  those  techniques  applied  to  the  scattering  from  rough 
solid  surfaces.  Application  of  these  ideas  to  wakes  was  first 
proposed  and  qualitatively  discussed  by  Salpeter  and  Treiman 
(Ref.  12).  The  analysis  of  the  electromagnetics  of  scattering 
from  rough  solid  surfaces  can  be  categorized  as  two  limiting 
classes  of  problems.  The  classes  are  called  either  slightly  or 
very  rough.  If  the  RMS  height  of  the  surface  roughness  is  much 
less  than  the  radar  wavelength,  the  surface  is  called  slightly 
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rough.  If  the  RMS  height  is  much  more  than  the  wavelength  it  is 
called  very  rough. 

For  slightly  rough  surfaces  one  usually  considers  the  scat¬ 
tering  cross  section  from  both  the  coherent  (specular)  and  in¬ 
coherent  (diffuse)  part.  The  coherent  part  is  that  which  comes 
from  a  smooth  surface.  The  incoherent  part  comes  from  the  rough¬ 
ness  itself.  Slightly  rough  surface  theory  is  derived  as  a  per¬ 
turbation  analysis  around  the  smooth  surface  case.  For  most 
problems  the  coherent  part  of  the  cross  section  turns  out  to  be 
much  greater  than  the  incoherent  part.  For  this  slightly  rough 
surface  coherent  cross  section  we  have  chosen  a  model  derived  by 
Peake  and  also  by  Davies  and  reported  in  the  radar  handbook  (Ref.  13). 
We  are  not  going  to  derive  the  slightly  rough  equation.  However, 
the  expressions  are  valid  under  the  condition  that 

kQh  <1  (18) 

where  h  is  the  RMS  roughness  height  and  kQ  is  the  wave  number. 

The  correction  factor  is 

Cr  =  exp  ~(VT  kQh  cos  6^)^  •  (19) 

The  expression  for  the.  incoherent  part  of  the  scattered  radiation 

is 

Cr  -  Akoh2/2  (20) 

where  the  coefficient  A  is  4  for  an  exponential  correlation 
function  and  8  for  a  gaussian  correlation  function,  and  /  is 
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the  surface  height  correlation  coefficient  of  the  surface  rough¬ 
ness  measured  along  the  surface.  We  will  see  later  that  the 
slightly  rough  electromagnetic  case  is  used  only  when  the  corre¬ 
lation  length  or  turbulence  scale  is  very  small.  For  this  case 
the  incoherent  scattered  radiation  is  much  less  than  the  coherent, 
and  therefore  it  will  be  neglected. 

For  increased  size  of  the  surface  roughness  compared  to  the 
radar  wavelength  one  finds  that  the  scattered  power  becomes  less 
and  less  coherent  until  at  some  large  RMS  roughness  height  the 
coherent  fraction  of  the  scattered  power  disappears  and  all  the 
scattered  power  is  then  incoherent.  The  surface  is  then  called 
very  rough.  The  very  rough  surface  correction  factors  are  valid 
under  the  conditions  that 

kQh  >  5  .  (21) 

The  assumptions  of  the  model  are  the  following:  a)  the  surface 
radius  of  curvature  must  be  larger  than  the  wavelength;  b)  iso¬ 
tropic  roughness;  c)  the  correlation  length  of  the  surface  rough¬ 
ness  must  be  less  than  the  scattering  area,  and  d)  multiple  scat¬ 
tering  and  shadowing  are  neglected.  For  backscacter  with  a 
gaussian  surface  correlation  function  the  roughness  correction 
factor  is 
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where  S  ■  4h  / l  and  is  an  average  angle  of  incidence. 

And  for  the  exponential  correlation  function 
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The  turbulent  fluid  scattering  surface  is  the  physical  area 
about  which  the  least  is  known  and  therefore  is  the  most  difficult 
to  model.  Two  different  physical  models  come  to  mind  (Pig.  13). 
Since  the  electron  density  decreases  radially,  it  is  possible  that 
the  surface  where  the  plasma  frequency  is  equal  to  radar  frequency 
could  be  interior  to  the  turbulent  plume.  The  region  interior  to 
the  scattering  surface  will  be  overdense.  This  model  of  the  scat¬ 
tering  surface  is  called  the  interior  turbulent  surface  model.  It 
is  assumed  that  the  surface  turbulent  properties  such  as  the  cor¬ 
relation  lengths  and  functions  would  be  the  same  as  the  volume 
properties  described  above.  Therefore  the  RMS  surface  height  and 
the  correlation  length  along  the  surface  are  equal  to  the  volume 
radial  and  longitudinal  correlation  lengths,  respectively.  Recent 
studies  (Refs.  5,  6,  7,  and  14)  have  shown  that  the  interior  tur¬ 
bulence  is  mainly  isotropic.  Therefore  the  correlation  length 
normal  to  and  on  the  surface  of  the  turbulence  will  be  equal.  For 


the  interior  turbulence,  calculations  were  made  for  plume  diameter 
to  correlation  length  ratios  (R/a)  of  10  (large  correlation 
length)  and  1000  (small  correlation  length) .  This  corresponds 


25 


to  correlation  length  of  0.97  and  0.0097  meters,  respectively. 

For  a  400  MHz  wave,  then,  from  expressions  (18)  and  (21)  the 
surfaces  are  slightly  rough  for  the  large  correlation  length  and 
very  rough  for  the  small  correlation  length.  For  the  interior 
turbulence  model  one  would  use  Eq.  (20)  for  the  slightly  rough 
roughness  correction  and  Eqs .  (22)  and  (23)  for  the  very  rough 
roughness  corrections. 

The  second  model  of  surface  turbulence  will  now  be  considered. 

In  many  cases  the  scattering  surface  is  the  surface  of  the  turbulent 
jet  itself.  This  jet  surface,  called  the  intermittent  front,  is 
undulating  in  a  statistical  manner.  The  intermittency ,  a  statistical 
parameter,  expresses  the  average  fraction  of  time  that  a  region  in 
space  near  the  front  is  interior  to  the  turbulent  plume.  There  are 
two  other  measures  of  the  intermittent  front  that  are  important  for 
electromagnetic  wave  scattering  (Ref.  15).  They  are  the  standard 
deviation  of  the  intermittent  front  radius  fluctuation  (o^), 
which  is  analogous  to  the  radial  correlation  length,  and  the  ap¬ 
parent  wavelength  of  the  turbulent  front  (A^)  ,  which  is  analogous 
the  axial  correlation  length. 

Demetriades  (Ref.  15)  has  measured  both  of  these  quantities  in 
a  plasma  jet  experiment  (Fig.  14).  His  jet  was  produced  by  passing 
argon  gas  through  a  dc  arc  and  expanding  it  through  a  supersonic 
nozzle  with  a  2 . 5  cm  exit  plane  diameter  into  a  large  vacuum  tank 
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held  at  10  torr  pressure.  Measurement  of  the  above  two  quanti¬ 
ties  of  interest  was  made  with  Langmuir  probes  and  hot  wire 
anemometers  in  the  turbulent  part  of  the  jet. 

Demetriades  reported  measurements  of  and  A  ^  with  re¬ 

spect  to  a  transverse  scale  length  L,  which  varies  linearly  with 
the  axial  position  of  the  turbulent  region  of  the  jet.  Furthermore, 
the  mean  diameter  of  the  turbulent  region  of  the  jet  increases  ap¬ 
proximately  linearly  with  increasing  axial  position.  The  value  of 
o^/L  is  1.5  and  the  values  of  A^/L  vary  between  a  maximum 
value  of  15  and  a  minimum  of  5.  From  this  one  can  calculate  a 
maximum  and  minimum  value  of  o/A^  =  h It.  of  0.3  and  0.1,  re¬ 
spectively.  Since  L  and  hence  both  of  and  A £  scale  linearly 
with  the  local  diameter  of  the  turbulent  portion  of  the  jet,  we 
assume  that  o^.Aph,  and  l  will  scale  linearly  to  the  full 
size  jet.  From  Eq.  (21)  we  find  that  all  the  intermittent  front 
cases  are  very  rough  surface  cases.  Therefore  the  surface  rough¬ 
ness  correction  factor  is  ted  from  Eqs.  (22)  and  (23). 

The  value  of  h il  =0.1  is  ca.  .«  large  axial  scale  length 

case  and  h/i  =  0.3  is  called  the  small  axial  scale  length  case. 

SURFACE  SCATTERING  FROM  A  SOLID  ROCKET  PLUME 


Eight  cross  sections  have  been  calculated  for  the  two  turbu¬ 
lent  surfaces  considered,  namely  interior  and  intermittent  front 


turbulent  surfaces.  All  results  are  backscatter,  broadside  cross 
sections.  The  interior  surfaces  cases  have  been  calculated  for 
the  large  correlation  length  with  both  gaussian  and  exponential 
correlation  functions  and  for  the  small  correlation  length.  All 
interior  surface  cross  sections  consider  isotropic  turbulence  only. 
For  the  interior  surface  case,  the  small  correlation  length  turbu¬ 
lence  leadi,  to  a  slightly  rough  electromagnetic  scattering  model 
and  the  large  correlation  length  to  a  very  rough  electromagnetic 
scattering  model. 

The  intermittent  front  cases  have  been  calculated  for  large 
radial  correlation  lengths  only.  However,  calculations  with  both 
large  and  small  axial- to-radial  correlation  length  ratios  and  both 
gaussian  and  exponential  correlation  functions  have  been  made.  All 
the  cases  treated  with  the  intermittent  front  models  lead  to  the 
very  rough  electromagnetic  scattering  models.  The  cross  section 
and  type  of  all  calculations  have  been  summarized  in  Fig.  15. 


CONCLUDING  REMARKS 

It  has  been  fully  recognized  for  many  years  that  both  the 
fluid  and  electromagnetic  aspects  of  the  study  of  rocket  plumes 
and  aircraft  jets  are  very  difficult  areas  of  engineering  physics. 
This  is  true  for  both  the  analytical  and  experimental  work  in  the 
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field.  In  this  paper  we  have  presented  the  calculated  radar  scat¬ 
tering  cross  sections  of  a  variety  of  analytical  models  of  both 
underdense  and  overdense  plumes.  Some  models  give  large  cross 
sections,  others  small  cross  sections.  The  models  include  physical 
phenomena  not  considered  in  detail  by  previous  investigators.  How¬ 
ever,  the  models  presented  by  no  means  exhaust  the  field.  One  can 
think  of  many  phenomena  that  might  bear  consideration  in  an 
analytical  model.  Two  things  can  be  noted,  however.  The  first  is 
that  the  "correct"  model  or  models  will  have  to  be  simple  in  order 
to  be  tractable.  Secondly,  the  models  must  also  explain  or  predict 
the  scattering  characteristics  of  plumes  more  accurately  and  in 
more  detail  than  older,  simpler  models  if  the  results  are  to  be 
useful  to  the  detection  and  discrimination  systems  of  the  future. 
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Boundary  of  Knee  in  Radial 
Electron  Density  Calculation 


y/R  *  10  from 
x  -  45  m  to  150 
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(Trail  Length  in  Ball  Diameters) 
Fig.  2  Correlation  Length  in  Wake  of  Sphere 
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Jet  Width/Correlation  Length 
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Fig.  4  One  Dimensional  Plasma  (Volume  Scattering) 


35 


2*2  +  xj  -  2A 


Fig.  5  Boundaries  of 
Space  for  the 


Fig.  6  Limits  of  Integration  for  I^ji  of  the 
Scattering  Cross  Section 


Fig.  8  Finite  Volume  to  Booker-Gordon  Scattering  Ratio 

versus  Plasma  Dimension  to  Correlation  Length  Ratio 
for  ka  =  0.1,  1,  2 


Correlation  Length  Ratio  for  ka 
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Fig.  11  Computer  Profile  of  Overdense  Region  for  Solid  Propellant  Plume 
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